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In this work we present a phase space analysis of a quintessence field and a perfect fluid trapped 
in a Randall-Sundrum's Braneworld of type 2. We consider a homogeneous but anisotropic Bianchi 
I (BI) brane geometry. Moreover, we consider the effect of the projection of the five-dimensional 
Weyl tensor onto the three-brane in the form of a negative Dark Radiation term. For the treatment 
£T") ■ of the potential we use the "Method of /-devisers" that allow investigating arbitrary potentials in a 

phase space. We present general conditions on the potential in order to obtain the stability of stan- 
dard 4D and non-standard 5D de Sitter solutions, and we provide the stability conditions for both 
■ scalar field-matter scaling solutions, scalar field-dark radiation solutions and scalar field-dominated 

solutions. We find that the shear-dominated solutions are unstable (particularly, contracting shear- 
dominated solutions are of saddle type). As a main difference with our previous work, the tradition- 
ally ever-expanding models could potentially re-collapse due to the negativity of the dark radiation. 
Additionally, our system admits a large class of static solutions that are of saddle type. This kind 
of solutions are important at intermediate stages in the evolution of the universe, since they allow 
to the transition from contracting to expanding models and vice versa. New features of our scenario 
are the existence of a bounce and a turnaround, which leads to cyclic behavior, that are not allowed 
in Bianchi I branes with positive dark radiation term. Finally, as specific examples we consider the 
a* potentials V oc sinh a (/3(f>) and V oc [cosh (£</>) — 1] which have simple /-devisers. 
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The idea that our Universe is a brane embedded in a 5-dimcnsional space, on which the particles of the standard 
model are confined, while gravity is allowed to propagate also in the bulk of the higher-dimensional manifold was 
, proposed in [3, 0. This braneworlds models were baptized as a Randall-Sundrum type 1 (RSI) and 2 (RS2) 1 . One 
of them, the RS2 model, found a great importance in cosmology and has been studied from different angles, because 
of the cosmological field equations on the brane have modification that have significant impact on inflationary 
scenario 0-0- 

^ ; The phase space of the Randall-Sundrum braneworld has been investigated widely in the literature El, Il3l4l7| . 
In [3 [l9| the asymptotic behavior of FRLW, Bianchi I and V brane in the presence of perfect fluid was studied 2 
taking into account the contribution of the no n- vanishing dark radiation termZY. Goheer and Duns by (20l.[2l| used the 
dynamical system methods to the analysis of how the FLRW and Bianchi I scalar field models with and exponential 
potential is affected by the bulk Weyl tensor. These studies was extended to arbitrary potentials in FRW branes 
j22l |23| while the contribution of the positive dark radiation term for arbitrary potentials in Bianchi I branes was 



INTRODUCTION 



generalized in [24 



The dark radiation term U represents the scalar component of the electric (Coulomb) part of the five-dimensional 
Weyl tensor of the bulk (2f|[26j . U scales just like radiation with a constant fj,, i.e., U = fia~ A . For that reason it is called 
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1 The original motivations were to looking for an explanation to the hierarchy problem Q|, in (RSI), and an alternative mechanism to 
the Kaluza-Klein compactifications Qj, in (RS2). 

2 A complementary analysis for FLRW branes was carry on in |10H using a 2-dimensional Hamiltonian system. 
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the dark radiation. The coefficient p is a constant of integration obtained by integrating the five-dimensional Einstein 
equations [1, 0, 127l-l32ll. In FRW brane model, p is related with the mass of the black hole in the bulk therefore U is 
positive defined [33l - !3a |. For anisotropic models both positive and negative p are possible mathematically fl9l.l2ll.l36l]. 
Its magnitude and sign can depend on the choice of initial conditions when solving the five-dimensional Einstein 
equation. Hence, even the sign of p remains an open question [261 [3~4| . 

Dark radiation should strongly affect both big-bang nucleosynthesis (BBN) and the cosmic microwave background 
(CMB). Hence, such observations can be used to constrain both the magnitude and sign of the dark radiation as well 
as any source of radiation in the model (see the classical references [37], [38[, see the appendix of [39[ and pol - l43l |). In 
[43j], the constraints on BBN alone allow for —1.23 < pdr/ Pj < 0.11 in a dark radiation component, where p 1 is the 
total energy density in background photons just before the BBN epoch at T = lMeV, pdr = 6K/X, where A is the 
tension of the brane. In order to compute the theoretical prediction of CMB anisotropics exactly, one must eventually 
solve the perturbations including the contribution from the bulk. However, in (43[ the CMB power spectrum was used 
to constrain the dominant expansion-rate effect of the dark-radiation term in the generalized Friedmann equation. If 
the constraint from effects of dark radiation on the expansion rate are included, the allowed concordance range range 
of dark radiation content reduces to —0.41 < pdr/ p-y < 0.105 at the 95% confidence level. So, in principle, the data 
allow for a negative extra radiation. Considering negative values of dark radiation in Bianchi I models leads to very 
interesting solutions for which the universe can both collapse or isotropize (l8l - [2ll |44| . All the above reasons motivate 
the investigation of a homogeneous but anisotropic brane model with U < 0. 

In this paper, we study the dynamics of an scalar field with arbitrary potential trapped in a Randall-Sundrum's 
Braneworld of type 2. We consider a homogeneous but anisotropic Bianchi I (BI) brane filled also with a perfect fluid. 
Furthermore, we consider the effect of the projection of the five-dimensional Weyl tensor onto the three-brane in the 



form of a negative Dark Radiation term. The results presented here complement our previous investigation in 24 1 
where we considered the effect of a positive Dark Radiation term on the brane. 

For the treatment of the potential we use a modification of the method that has been applied to isotropic (FRW) 
scenarios [22|, I45T - I491 ] , and that has been generalized to several cosmological contexts [H, [2J, [5(1 HJ . The modified 
method, that we call "Method of /-devisers" , allows us to perform a phase-space analysis of a cosmological model, 
without the need for specifying the potential. This is a significant advantage, since one can first perform the analysis 
for arbitrary potentials and then just substitute the desired forms, instead of repeating the whole procedure for every 
distinct potential. This investigation represents a further step in a series of works devoted to the use of the gen eral 
procedure of /-devisers for investigating both FRW and Bianchi I brancs, initiated in our previous works [23L |24|. 

A main difference with our results in [24j , is that the traditionally ever-expanding models could potentially rccollapse 
due to the negativity of the dark radiation (U < 0). Also, our system admits a large class of static solutions that 
are of saddle type. This kind of solutions are important at intermediated stages in the evolution of the universe since 
they allow for the transition from expanding to contracting models, and vice versa. New features of our scenario are 
the possibility of a bounce and a turnaround, which leads to cyclic behavior. This behavior is not allowed in Bianchi 
I branes with positive dark radiation term [24j . 

The interest of static solutions in the cosmological setting goes back to the discovering of Einstein static (ES) 
model. This model was proposed by Einstein as an attempt to incorporate Mach's principle (geometry or motion do 
not make sense in an empty universe) into the general relativity (GR) and also to overcome the boundary conditions 
of the theory 
by Eddington 
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The issue of the stability of ES model has been studied several times after the classical paper 
53| where it was showed the instability of ES solution with respect to the homogeneous and isotropic 
perturbations. This fact is the key feature that allows for the transition between the Friedmann transient matter- 
dominated deceleration phase, to the current accelerated expansion. In [54| . and [55| . were considered inhomogeneous 
and anisotropic perturbations to ES model. A cosmology including a mixture of isotropic radiation and a ghost 
masslcss scalar field admits an almost static solution that it is stable against homogeneous an isotropic perturbations. 
These solutions are also stable against the three types of linear inhomogeneous perturbations: density, vorticity 
and gravitational-wave perturbations. However, ES model becomes unstable (it no longer exists) if large-amplitude 
homogeneous and anisotropic perturbations are introduced [56j . The stability issue of ES models has acquired renewed 
interest, for example, in the context of f{R) theories due the success of such theories in explaining the cosmic 
acceleration and because ES solutions can provide the link between decelerating/accelerating phases in these theories 
in the same way as in GR [13] ■ The stability analysis in [57j shows that as in GR [58|, the Einstein static model is 
neutrally stable against vector and tensor perturbations for all equations of state at all scales. In (59[, is presented 
a complete picture of how ES models fits into the overall phase space structure of f{R) gravity, complementing the 
perturbation analysis made in their previous wo rk [57l| . There the authors investigate the prototype of /(i?)-theory 
compatible with an static solution obtained in [57] . Several unexpected cosmological features were obtained such 
as the existence of no big-bang scenario for any realistic equation of state, but only bounce scenarios or expansion 
after an initial Minkowski phase; the existence of past and a future bounce, a phenomenon which is absent in general 
relativity; the unfolding of the quasi-periodic nature of ES model, in the sense that orbits emerging from a source, 
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spiral around ES solution for some time and goes to a sink for a negative equation of state parameter of the perfect 
fluid. The stability analysis of the ES universe, and other kinds of static solutions in the context of brane cosmology 
with dust matter and anisotropic geometry was presented in [lj| . Here we present several classes of static solution 
for a Bianchi I brane containing an scalar field with arbitrary potential. This class of solutions are more general than 
the presented in in Il9j . As commented before, our static solutions and very important since they allow to the 
transition from a contracting to an expanding universe, and vice versa. 

The paper is organized as follows. In section [TT] we present our general approach for the investigation of arbitrary 
potentials. In section IIIII are presented the cosmological equations of our model. In section IIVI we proceed to the 
dynamical system analysis of the cosmological model under consideration. For this purpose we use the method of 
/- devisers discussed in section [III an d we introduce local charts adapted to each of the more interesting singular 
points. In this way, we analyze the stability of the shear-dominated solutions, the stability of the de Sitter solutions 
and the stability of the solutions with 5D corrections. In section fVl we summarize our analytical results and proceed 
to the physical discussion specially to the comparison with previous results. Particular emphasis is made on static 
solutions. Our analytical results are illustrated in section |Vl] for the potentials V(6) = Vb [cosh (£^>) — 1] (see references 
H3, SB, IH, ) and V((j>) = Vosinh~ a ((3<f>) (see references 0, H ED, E3-EE Ezl ) which have simple /-devisers. 

Is worthy to mention that this procedure is general and applies to other potentials different from the cosh-like, the 
inverted sinh-like and the exponential one. Finally, in section fVIII are drawn our general conclusions. 



II. "METHOD OF /-DEVISERS 



In this section we present a general method that allows us to perform a phase-space analysis of a cosmological 
model, without the need for specifying the potential. This procedure that we will call "Method of /-devisers" consist 
of a modification of a method used in isotropic (FRW) scalar field cosmologies for the treatment of the potential 
of the scalar field [22], I45l - l48j , that has been generalized to several cosmological contexts 0, [24], [5(1 [H| ■ An earlier 
attempt to introduce the "Method of /-devisers" as in the present form was in (4(| in the context of FRW cosmologies. 
However, there the authors restricted their attention to a scalar field dark-matter model with a cosh-like potential, 
and there was revealed that the late-time attractor is always the de Sitter solution. In section I VI Al we will take 
this potential as an application example of the present "Method of /-devisers" in the context of brane cosmology. 
The method as introduced here applies to a more general class of potentials than the cosh-like. In fact, the method 
has the significant advantage that one can first perform the analysis for arbitrary potentials and then just substitute 
the desired forms, instead of repeating the whole procedure for every given potential. More importantly is that the 
method does not depend on the cosmological scenario, since it is constructed on the scalar field and its self-interacting 
potential, and can be generalized to several scalar fields (for example to quintom models [H| ) . The disadvantages are 
that such a generalization is only possible if in the model do not appear functions that contains mixed terms of several 
of the scalar fields. For example, cosmological models containing an interaction term given by V((j), ip) discussed in 
references [6^ | and [7^, where <fi,tp are the scalar fields, cannot be investigated using the method of /-devisers. On 
the other hand, in the single field case the method cannot be applied if the model contains more than one "arbitrary" 
function depending on the scalar field, as for example a cosmological models containing a potential V(4>) and coupling 
function x(<^)i that appeared in the context of non-minimally coupled scalar field cosmologies [7ll - l73l ]. In such a case 
the more convenient approach is to consider the scalar field itself as a dynamical variable. 

The procedure for deriving the "Method of /-devisers" is as follows. Let us define the two dynamical variables 



V{4>) 
V{4>) 



(i) 



1 v(<t>) V{W ' 1 > 

while keeping the potential still arbitrary. 

In the formulas ((T|) and ([2]), the primes denote differentiation with respect to (j>. Now, a simple observation is 
that for general physical potentials, and for the usual ansatzes of the cosmological literature, it results that / can be 
expressed as an explicit one-valued function of s, that is / = f(s) (as can be seen in Tablc[T|, and therefore we result 
to a closed dynamical system for s and a set of normalizcd-variablcs. 

Conversely, it is possible to reconstruct the corresponding potential from a given f(s). In particular, one starts 
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TABLE I: Explicit for the function f(s) for some quintessence potentials. 



Potential 


References 


/(«) 


V(<p) = V [cosh^<j>)-l] 


[22, 46, 48, 60-66] 




V(<f>) = V smh- a (l3<t>) 


[22, 48, 61, 63-65, 67] 




V(<t>) = V e- x<t, + A 


[63, 74J 


-s(s-\) 


V(4>) = V [e Q< * + e^] 


[75-77] 


-(s + a)(s + f3) 



with 



£ = -V. (4) 



which give 



<Ks) = <h- I jm dK ( 5 ) 



so 



V(s) = e J ^7Wl dK V , (6) 

where the integration constants satisfies V(sq) — Vq, 4>(so) = 4>a- We mention that relations ((5]) and (j6|) are always 
valid and they provide the potential in an implicit form. However, for the usual cosmological cases of Table U the 
potential can be written explicitly, that is V = V (</)), after elimination of s between ([5]) and (|6]). 

The direct derivation of the function f(s) from a given V((j>), as well as the reconstruction procedure for obtaining 
V(4>) from a given /-function is what we call "Method of /-devisers" . In the section IIVI we use this method for 
investigating the phase portrait of quintessence fields trapped in a Randall-Sundrum Brancworld in the case of an 
anisotropic Bianchi I brane with a negative dark radiation term. 



III. BIANCHI COSMOLOGY ON THE BRANE AND QUALITATIVE BEHAVIOR 

Having presented the method we will use in the present paper, now let us present the cosmological model. We will 
be concerned here with the dynamics of a self interacting scalar field with an arbitrary self-interaction potential, that 
is trapped in a Randall-Sundrum brane of type 2 (RS2). The field equations for the Bianchi I metric are the following: 

3 



H 



3 PT [ 1+ 2X) + ~f 



•1U 



(7) 



* = -i(l + f)(* +m 



a = -3Ha 



4U 



p rn + 3H(p m +p m ) =0 



(8) 

(9) 
(10) 



3 See our previous paper [531 and and therein references for a full derivations of the effective Einstein equations for Bianchi I models. 
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^ + 3fty + 00V = O, (11) 

where we have used units in which k 2 = 8(f>G = 1. In the above expressions denotes the Hubble factor, a is the 
shear term, while 

PT = P<f> + Pm (12) 

is the total energy density content of the universe, p m corresponds to the energy density of a barotropic fluid with 
equation of state 

Pm = (7-l)Pm, (13) 

and 

p^ y + V(<t>) (14) 

is the energy density of the scalar field <j>. U denotes the so-called Dark Radiation term which in this paper is assumed 
to be negative, Li < 0, as it could be the case according to observations [37| H3|. From the equations and (|5J), and 
from the fact that U < 0, it follows that H can become zero and moreover change its sign during the evolution. Thus, 
for U < the expanding phase could experience the cosmological turnaround triggered by the negativity of this dark 
radiation, and similarly the contracting phase can be led to a cosmological bounce and then to expansion. 

Before proceeding to the detailed investigation of the model (f7])- (fTTj) using dynamical systems tools, let us discuss 
the above system using a heuristic reasoning. Integrating equation ([9]) and ()10[) we obtain that a = troa 3 and 

Pm = PmO a ~ 37 j while as discussed in the introduction U = fj,a~ 4 . Let us examine the case of 4- S> V (</)), thus, 
P<p ~ P0oa _6 i which implies pt = p^oa^ 6 + /9 m oa~ 37 . Then equation ([7]) reduces to 

Qu2_ PmO a ~ 67 | PmOP0OQ~ 37 ~ 6 flfl) -12 , /. , _2\ 6 , _ „-3 7 6 ImI 4 , K1 

d-H — 1 V — a +(P0O + cr oJ a +Pmoa —a (15) 

There are several cases depending of the values of 7. Suppose that the universe lies currently in the usual expanding 
phase, i.e. H(a = ao) > 0, where ao denotes the scale of the universe at the present time. We take as a reference 
universe state the present expanding one. Thus a future zero of H will correspond to a turnaround, while a past one 
corresponds to a bounce. 



For 7 > 4/3 the dark radiation term in equation (|15|) falls less rapidly than the matter terms, making the future 
turnaround inevitable for any initial negative 14. 



• For 7 = 4/3, the last two terms in the r.h.s. of (|15[) are equally important in a low-energy regime, and the 
condition for the future turnaround is p m o < 6\p\/X. 

• For 1 < 7 < |, the leading terms in the r.h.s. of (|15|) that are equally important are the last two terms (since 

the terms scaling as a -6 goes very fast to zero). As a result, we obtain that the condition H = leads to the 

1 

equation /9 m oa 4 ~ 37 — = 0. This solution as a positive root given by a+ = T . If a+ > ao (the 

scale factor today) it corresponds to the future of the Universe under consideration and indicates the point of 
turnaround. If a+ < ao, will describe the past and corresponds to a bounce. 

• For 7=1, the leading terms in the r.h.s. of (|15JI that are equally important are the first and the last three 
terms. As a result, we obtain that the condition H = leads to the cubic equation for the scale factor 

fl3 ~ Xp^o 0,2 p^o ""^o = ^' APPh/ m § Descartes's rule, this equation has always one negative root 

and two real positive or either complex conjugated roots. The condition for the existence of two positive roots 
is 



PmO + 2 Vmo(P0O + CTo) < ' 

This expression reduces to the analogous expression in Q for the case of a perfect fluid under the assumption 
P4>o = <7o = 0. A root which is bigger than ao (the scale factor today) corresponds to the future of the Universe 
under consideration and indicates the point of turnaround. A root which is less than ao, will describe the past 
and corresponds to a bounce. 
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• For | < 7 < 1, the leading terms in the r.h.s. of (|15|) that are equally important are the first and the last two 
terms. As a result, we obtain that the condition H = leads to the equation 

t r p,noa + - — U. 

This equation must be solved numerically looking for real positive roots. A root which is bigger than ao (the 
scale factor today) corresponds to the future of the Universe under consideration and indicates the point of 
turnaround. A root which is less than ao, will describe the past and corresponds to a bounce. 

• For 7 = | , the leading terms in the r.h.s. of (1151) that are equally important arc the first and the last two 

leads to the quadratic equation for the scale factor 

, is PmO < 12 N- 



terms. As a result, we obtain that the condition H 

' (y^TT ~ Xp^o ) = ^' coll dition for the existence of a positive root, a + = (^y 1 



a 



6H 

Ap m C 



If a + > ao (the scale factor today) corresponds to the future of the Universe under consideration and indicates 
the point of turnaround. If a + < ao, it describes the past and corresponds to a bounce. 

• For 7 < 2/3, even the term oc a -67 , falls less rapidly than U, and the turnaround in the future becomes 
impossible. 



Let us examine the case of ^- -C V (</>), thus, 
reduces to 



3H Z = V 



1 



V_ 

2A 



PmO a 

2A 



V{4>{a)). Then, p T = V{4>{a)) + p m0 a -37 ■ Then equation © 

(16) 



PmO 



1 



V 

7 



,-3 7 



6N 
a 4 A' 



Let consider the simpler case of a potential that tends asymptotically to a positive constant V = Vq > 0. 

• For 7 > |, the leading terms in the r.h.s. of ([TS]) are the first and the last one. The condition H = leads to 



the quartic equation for the scale factor 



V 



6H 
A 



This equation has one negative solution, two purely imaginary solutions, and a positive one given by a+ = 

\P1 v^3 4 /ZT 

, v . If a + > ao (the scale factor today), it corresponds to the future of the Universe under consideration 

V Vb(2A+Vb) 

and indicates the point of turnaround. If a+ < ao, it describes the past and corresponds to a bounce. 

For 7 = f j the leading terms in the r.h.s. of (Tip]) are the first and the last two terms. The condition H = 
leads to the quartic equation for the scale factor 



Va 



PmO 



1+^ 



6M 



= 0. 



The condition for the future collapse is p m Q [l + ^-1 < 5M, 



• For 1 < 7 < |, the leading terms in the r.h.s. of (JT6J) are the first and the last two terms. The condition H 
leads to the equation for the scale factor 

Yi 

A 



V 



1 + 



2A 



PmO 



1 



,4-37 



6_H 
A 



0. 



This equation must be solved numerically looking for real positive roots. A root which is bigger than ao (the 
scale factor today) corresponds to the future of the Universe under consideration and indicates the point of 
turnaround. A root which is less than ao, will describe the past and corresponds to a bounce. 

• For 7 = 1, the leading terms in the r.h.s. of (fT6|) are the first and the last two terms. The condition H = 
leads to the quartic equation for the scale factor 



Vo 



PmO 



6H 
A 



= 0. 



Applying Descartes's rule this equation has one positive root, one negative root and two complex conjugated 
roots. If the positive root is bigger than ao (the scale factor today), it corresponds to the future of the Universe 
under consideration and indicates the point of turnaround. If it is less than ao, it describes the past and 
corresponds to a bounce. 
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For | < 7 < 1, all the terms in the r.h.s. of (|T6|) but the third, are relevant. The condition H = leads to the 
equation for the scale factor 



6/xa 6 ^ 4 , p m0 a^{X + V ) , V a^(2X + V ) 



X 



X 



2X 



2 

PmV) 

2X 



0. 



This equation must be solved numerically looking for real positive roots. A root which is bigger than ao (the 
scale factor today) corresponds to the future of the Universe under consideration and indicates the point of 
turnaround. A root which is less than ao, will describe the past and corresponds to a bounce. 

For 7 = |, all the terms in the r.h.s. of (fT5| but the third, are relevant. The condition H = leads to the 
quartic equation for the scale factor 



Vo 



1 



Vb 
2A 



PmO 



X 



2 

PmO 

2X 



6\fA 
X 



= 0. 



For [i > ^yrp there is one positive solution given by a + 



V x2 P 2 m + Vo (2 A+ V a ) +p m0 (- A- V ) 
V (2A+Vb) 



, one negative 



solution a_ = — 



V A2 PmO + 12 A'Vc)(2A+Vo)+p m o(-A-yo) j , , , , • r 

x — Vo(2A+y ) ' a com P^ ex conjugated. Otherwise the four are 



complex conjugated. Thus, the condition for the bounce or the collapse is fi > ^-ff-- If the positive root is bigger 
than ao (the scale factor today), it corresponds to the future of the Universe under consideration and indicates 
the point of turnaround. If it is less than ao, it describes the past and corresponds to a bounce. 

• For 7 < 2/3, even the term oc a -67 , falls less rapidly than U, and the turnaround in the future becomes 
impossible. 

Summarizing, using a heuristic reasoning we have obtained several conditions for the existence of bounces and 
turnaround of cosmological solutions for a massless scalar field and a scalar field with a potential asymptotically 
constant. For the analysis of more general potentials we submit the reader to the reference fll| where the authors 
discuss conditions for bounces and they provide details of chaotic behavior in the isotropic brane model containing 
an scalar field for a negative dark radiation U < 0. 



IV. DETAILED DYNAMICAL ANALYSIS FOR BIANCHI I BRANE WITH NEGATIVE DARK 

RADIATION TERM 

After the qualitative analysis of the previous section, in the present section we investigate the cosmological model 
(|7])- ([TTj) from the dynamical systems perspective. First we need to recast the cosmological equations (|7|-(fTT|) as an 
autonomous system of first order ordinary differential equations (72l . l73l frsMsol ] . In order to avoid ambiguities related 
to the non-compactness at infinity we define compact variables that can describe both expanding and collapsing 
models [H GJ, ELr® : 



H V 



Pt 



D' 3£P' 6XD 2 ' 

°— ^• E -VSD' lfc — ^ ,17 » 



where D = J H 2 - 

For the scalar potential treatment, we proceed following the method introduced in section HTl 
Using the Friedmann equation we obtain the following relation between the variables (JT7J) 

x 2 + y + a m + n x + s 2 = i (is) 

The restriction (|18p allows to forget about one of the dynamical variables, e.g., y, to obtaining a reduced dynamical 
system. From the condition < y < 1 we have the following inequality 

< x 2 + Q, m + n x + S 2 < 1 (19) 
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The definition of D leads to the additional restriction 

Q 2 + n lr = l. (20) 



Thus, restriction (|20j) allows to eliminate another degree of freedom, namely variable flu- 

Using the variables (|17[) . the field equations (|71)- (fTT|) and the new time variable dr — Ddt, we obtain the following 
autonomous system of ordinary differential equations (ODE): 

Q'=\ (Q 2 - 1) (3 T Q TO + 6E 2 + 6x 2 - 4) + ^3 (Q 2 - l) (21) 

x' = i (?>Qx ( 7 O m + 2£ 2 + 2a- 2 - 2) + \/6s (l - x 2 - fl m - fi A - £ 2 )) + ^EQx (22) 

0^ = 3QO m ( 7 (0 m - 1) + 2£ 2 + 2x 2 ) + 3SQO m (23) 
fl' x = 3Qil x ( 7 Q m + 2S 2 + 2.x 2 ) + 3SQ(ft A - 1) (24) 

£' = ^QE ( 7 fi m + 2S 2 + 2.x 2 - 2) + ^SQ£ (25) 

s' = -VQxf(s). (26) 

Where the comma denotes derivatives with respect to r, and 

S = ^ ( 7 O m + 2a; 2 ) , (27) 

From (fT5|) follows the useful relationship 

PT 2ft A 2f* A 



A x 2 +y + fl m l-f7 A -£ 2 
It is easy to sec from (|2"5)l and the definition ([2T| that 

„ _ 2fl x ( 7 » m + 2x 2 ) 
" _ 1 - fix - £ 2 



(28) 



(29) 



From (|28p follows that the region 1 — A — £ 2 = TO +a; 2 +j/ = corresponds to cosmological solutions where pt S> A 
(corresponding to the formal limit A — > 0). Therefore, they are associate to high energy regions, i.e., to cosmological 
solutions in a neighborhood of the initial singularity 4 . Due to its classic nature, our model is not appropriate to 
describing the dynamics near the initial singularity, where quantum effects appear. However, from the mathematical 
viewpoint, these region + S 2 = 1) is reached asymptotically. In fact, as some numerical integrations corroborate, 
there exists an open set of orbits in the phase interior that tends to the boundary f2 A + £ 2 = 1 as r — > — oo. Therefore, 
for mathematical motivations it is common to attach the boundary il\ + E 2 = 1 to the phase space 5 . On the other 
hand the points with (Q A = 0) arc associated to the standard 4D behavior (pr <C A or A — > oo) and corresponds to 
the low energy regime. 

From definition (|17[) and from the restriction (|18p. and taking into account the previous statements, it is enough 
to investigate to the flow of (|2"2"|) - (|2"6"|) defined in the phase space 

* = {{Q,x,n m ,n x ,Z) : < a; 2 + ft m + fl x + E 2 < 1,-1 < Q < 1, 

-1 < x < 1,0 < O m < 1,0 < Q A < 1,-1 < S < 1} x {s G M}. (30) 

The system (|2"2")) - (f2l)]) admits twenty four classes of (curves of) fixed points associated to expanding (contracting 
solutions). Its coordinates in the phase space are given in Table HI Note that for U = are recovered the models that 
we have been investigated previously in (23|. As we commented previously, our model is not applicable near initial 
singularity (that is at the singular surface 1 — fl x — E 2 =0). However, in principle we can apply a similar approach 
as in [2~H to investigate the dynamics on the singular surface. We submit the reader to the Appendix [B] for further 
details. 



4 See the references [7ll. l87t for a classical treatment of cosmological solutions near the initial singularity in FRW cosmologies. 

5 We submit the reader to the Appendix [E] for a discussion of how to dealing with the singular points at the hypcrsurface £l\ + E 2 = 1. 



9 



TABLE II: Critical points of the system (I21I) -(I26 [ ) associated to expanding (contracting) solutions and their existence conditions. 
We use the notation s* for an s- value such that f(s*) = and s c for an arbitrary s- value. 



Label 



Qt 

Q±( S *) 

Qi(**) 
Qi(s*) 

Qt 

Qs(«*) 
Qt 
Qio 
Qfi 

Q«(0 



Q 


.(' 




tlx 


£ 


s 


Existence 


±1 





l 








s c g R 


always 


±1 


1 











s* 


always 


±i 


-1 











s* 


always 


±1 


y 2 s* 


I _ 3 T 








s* 


s* 2 > 37 


±i 













s* 


s* 2 < 6 


±i 











-1 


s c G R 


always 


±i 











1 


s c G R 


always 


±i 


cos u 








sin-u 


s* 


always 


±1 

















always 


±i 








«Ac G (0,1) 








always 


±1 








1 





s c G R 


always 


- 1 - 2 













s* 


s* 2 < 4 



The eigenvalues of the linear perturbation matrix evaluated at each of these critical points are displayed in the 
table HIl 

In section HV Al are discussed the stability conditions for the expanding (contracting) solutions. 

Additionally the system (|22T) - ([26l) admits eighteen classes of (curves of) fixed points corresponding to static solutions 
(H = Q = 0). They are displayed in table HVl The eigenvalues of the linear perturbation matrix evaluated at each of 
these critical points are displayed in the table [V] In section IIVBI are discussed the stability conditions and physical 
interpretation for this static solutions. 

Summarizing the system (pZl?f - (P()f admits forty two classes of fixed points. For that reason this scenario has a 
rich cosmological behavior, including the transition from contracting to expanding solutions and vice versa. Also is 
possible the existence of bouncing solutions and a turnaround, and even cyclic behavior. 



A. Stability of expanding (contracting) solutions 

Now let us comment on the stability of the first order perturbations of (|21 |) - ([26|) near the critical points showed in 
table El 

The line of singular points , although it is non- hyperbolic (actually, normally hyperbolic), they behave like a 
saddle point in the phase space of the RS model, since they have both nonempty stable and unstable manifolds (see 
the table [TTTj ) . The class Qi~ is the analogous to the line denoted by P\ in [24| . 

The singular points Q^(s*) and Q^(s*) are non- hyperbolic (these points are related to investigated in [HI), 
however they behave as saddle points since they have both nonempty stable and unstable manifolds (see the table 

ED. 

The singular point Q^(s*) is the analogous to P4 in (24[. It is a stable node in the cases < 7 < |,s* < 
-V5y,f (O < or § < 7 < 3,-7^5 ^ s * < -V^f,f(s*) < 0, or < 7 < §,** > V^f, f («*) > 0, or 
|<7<!,^< S *<^,/'( S *)>0. 

It is a spiral stable point for cither § < 7 < §, s* < --^t^J' (s*) < or § < 7 < |, s* > -j=^,f 0*) > 0. 

In summary, it is stable for either < 7 < |, s* < —^/3j, f (s*) < or < 7 < |, s* > y/Wj, /' (s*) > 0. Otherwise, 
it is a saddle point. 

The singular point Q4 ( s *) is a local source under the same conditions for which Q^(s*) is stable. 
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TABLE III: Eigenvalues for the critical points in Table [TTJ We use the notations (7,a*) 



_3 ^ _ 7 ± V^^H2-9,)+2^ _ L±(fiA) = _i (3 ± ^12/(0)^-1) +9) and M±(a*) = 1 (-a* ± ^64 - 15 (a*) 2 
where s* denotes an s-value such that f(s*) = 0. 

Label Eigenvalues 

Qt {(3 7 - 4)c, |( 7 - 2)c, |(7 - 2)e, -3 7 e, 3 7 6, 0} 

Qa («*) {2e, 0, -6e, -3( 7 - 2)e, 6e - y/§8*,-y/§f (a*)} 

Q±(a*) {2e, 0, -6e, -3( 7 - 2)e, 66^/fe*, >/6/' (a*)} 

Qf (a*) {-37^, |(7 - 2)6, (3 7 - 4)e, eK+ (7, a*), eK~ (7, O. -^^} 

{-6 ( S *) 2 , 1 (a* 2 - 6) e, § (a* 2 - 6) 6, (a* 2 - 3 7 ) 6, (a* 2 - 4)6, -6a*/' (a*)} 
Q± {66,66,26,0,0,3(2 -7)6} 

Q7 {66,66,26,0,0,3(2 -7)6} 

Qf(s*) {0, 6e - v^a cos(w), -3(7 - 2)e, 2e, -6e, — v/6 cos(w)/'(s)} 

Q9 {-46, -36, 0, 0, 1 (-v/9 - 12/(0)3) 6, 1 (y9 - 12/(0) - 3) 

{-46, -36, 0, -3 7 6fi A , 6L+(f2 A ), 6L-(fi A )} 
Qn {2(3 7 - 2)6, 3(7 - 1)6, 3 7 6, 3 7 6, 3( 7 - 1)6, 0} 

Q±(a*) {-2 S *6,-V,i(4-3 7 )a*6,6M+( S *),6M-(a*),-26/'( S *)} 



The singular point Q§(s*) is the analogous to P5 in [24]] . It is not hyperbolic for s* € {0, ±\/6, ±\/37, 2} or 
/'(«*) =0. 

In the hyperbolic case, Q^s*) is a stable node for either < 7 < |, < s* < 0, /' (s*) < or | < 7 < 2, — 2 < 

s* < 0, f (s*) < or < 7 < f , < s* < V37, /' (s*) > or § < 7 < 2, < s* < 2, /' (s*) > 0; otherwise, it is a 
saddle point. 

The singular point Q$(s*) is a local source under the same conditions for which Qg~(s*) is stable. 

For evaluating the Jacobian matrix, and obtaining their eigenvalues at the singular points Qq, Q^ and we 
need to take the appropriate limits. The details on the stability analysis for Q§ and are offered in the Appendix 
(|B|) . The analysis of Qy is essentially the same as for Qg . 

The circles of critical points Q s are non- hyperbolic. Both solutions represents transient states in the evolution of 
the universe. 

The singular point Qg and the line of critical points Qtn are non-hyperbolic (Qg is the analogous of P±o for 
fix = while Qf is the analogous of P10 for 0> 7^ in [24[). In the Appendix lAl we analyze the stability of Qg by 
introducing a local set of coordinates adapted to the singular points. The non-hyperbolic fixed point Qg , behaves as 
a past attractor. 

The singular point Q^ is the analogous of Pu in [24|, this point represents a ID set of singular points such that 

= 1), parameterized by the values of s c . It is normally hyperbolic. 
The singular points Q^ 2 (s*) are non-hyperbolic for s* = 0,±2, 7 = 4/3 or f'(s*) = 0. Qi 2 (s*) is a stable node 
f<7<2,2<s*< ^7=,/' (s*) > 0, a stable spiral for § < 7 < 2,s* > -^=,/'(s*) > 0. In summary, Qf 2 (s*) is 

stable for | < 7 < 2,s* > 2,/' (s*) > 0. Otherwise it is a saddle point. <2^ 2 ( s *) * s unstable for the same conditions 
for which Qf 2 (s*) is stable. 

Summarizing, the possible late time attractors are: 

• Qt(s*) is stable for either < 7 < f , s* < -y/frj, f (s*) < or < 7 < f , s* > ^37, /' (s*) > 0; 

• Qf(s*) is stable for cither < 7 < f,->/37 < s* < 0,/'(s*) < or | < 7 < 2, -2 < s* < 0, f (s*) < or 
0<7<|,0<s*< V3^, /' (s*) > or § < 7 < 2, < s* < 2, /' («*) > 0; 

• and the line Q10 are stable for /(0) > 0. However, if we include the variable y = 3^7, in the analysis, the 
solutions associated to the line of fixed points Q^ are unstable to perturbations along the ?/-direction. 



TABLE IV: Critical points of the system (121[) - (126[) representing static solutions and their existence conditions. We use the notation s* for an s-value such that f(s*) = 

2^3S^ — 2^ ( + — 1^ a fjv^ o — 2S^ — — 2S^ 

and s c for an arbitrary s-value. We use the notation fi(-y, Q A , E) = 37 (o A _s 2 +i) — ' ^ 2 (^ A ' S ) = =3£*+3sT^+5 > / 3 (t> S) = 2 (l - E ) - ^— - — , / 4 (7, E) = ^— h 



E 2 -l, / 6 (E) = ff^f, / 6 ( 7 ,z,E) = 



-6E 2 +6a; 2 -y4(3(7-l)S 2 -37+2) :4 +9(7-2) 2 i 4 + 12(7-2)i 2 (3(7+l)S 2 -37-2)+4 6z 2 + V4(2-37) 2 +9(7-2) 2 3 : 4 + 12((4-37)7+4) : c" 2 "+4 

67 ' J7l7j a; J — 67 



Label Q 


X 




n x 


s 


s 


Existence 



Ef n Ac e [0, |] ±J§ < 7 < 2 

£2 cos 2 u sinw n Ac fg {±1, ± ^/| } , < 7 < 2, < ti < 2ir 

Ea /i(7,n Ac ,S c ) n Ac e[0,l-E 2 ) £ c 6 (-\/f, \/f) f 2 (n Xc , E c ) < 7 < 2 

£4 /s(7,E c ) /4(7,E C ) E c s* -l<Ec<-\/f,0<7<2or 



3 <E C < Jf,/ 6 (Sc)<7<2or 



< £ c < 1,0 < 7 < 2. 



x c l-E 2 -x 2 -ft Ac -\ +/ 6 (7,x c ,E c ) E c s* 2 2 + E 2 < §, < 7 < 2 



2 



2 



Ef ±«/| - £ 2 E c — t/f <E C < v/f,0<7<2 

S± 1-2E 2 E c E 2 < |,E C ^ 0,0 < 7 < 2 

B 8 x c l-x 2 c -fl Xc _^ + / 7 ( 7jXc ) s* x\ < |,0 <7< i^gl 

£9 s^o^+i) ^ c < 7 < §,n Ac = lor 

. <7< <«Ao<l, or 

I < 7 < 2,0 < n Ac < 1. 

4 3lS) !- 2S c S c s* 7^|,E c = 0or 

< 7 < |,0 < E 2 < \. 

-En 1-2E 2 E c 0<7<|,E c =0or 

f <7<2,E 2 <1. 

Bia 2(2 3 7 Sg) E c 0< 7 <|,l-3(2 2 3^<^<|or|<7<2,E 2 <|. 



±V3& 3(2^) S * 0< 7 <|. 
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TABLE V: Eigenvalues for the critical points in Table IIVI We use the notations 51(7, fix,E) = 

v/-3 7 (3E2-2)(-E2 + f2 A + l) 2 + 18E6-44S 4 +2E 2 (n A (9n A +2) + 17)-S(n2 + i) \/7(0) V s2 + n A "1 v /3 7 (-S 2 +n x + l)+6S 2 -4 

_ &+Qx+1 , fla(7,nx,E) ^ V , E2+nA+1 

93(7, x, fix, E) -§ (7 2 (fiA - l)fi. + (7 - 2) 2 S 4 - (7 - 2)S 2 (2 7 fix + 7 - 2)) D + n 2 , + i) 

9(7-2)a: 2 (-0 A (4(7-l)E 2 +7+2) + (7-2)(2E 4 -3E 2 + l)+27n 2 - ) 
2(E 2 +fl A -l) ' 
3j-7 2 (nx-l) 3 ^-(7-2) 2 x"l(n 2 +l)-(7-2) a: 2 (r! > -l)(27n 2 ,-(7+2)n A +7-2) 

94(7, fix) = 72 ( n A -i) 

n , s _ \/67(n A +i) 2 -8(n 2 +i) _ v 7co) A /fwy3 7 (n A +i)-4 
5s( 7 , fix) - fj^pj , 36(7, «a) - V ^n. + i ' 

97(7, E) = %/V 7 (9E 2 - 6) - 4E 2 + 4, fl8 ( 7 , E) = .Ay (6 - 9E 2 ) + 18E 2 - 8, 



and. 99 (7,E)= ^ 3( - 2)S2 - 3 - +4 . 

Ai(e, E c ), A2(e, E c ), A 3 (e, E c ) and A4(e, E c ) are the roots of the equation 
A 4 + AV4-6E 2 e/'(0) + A 2 [/(0) - 4] - 4AV4 - 6Ege/'(0) - 6/(0)E 2 = 0. 

fj, 1 (e,E c ), f i 2 (e, E c ),/i 3 (e, E c ) and M4(e,S c ) are the roots of the equation /i 4 + ^2^ 3 E c e/'(0) + fi 2 [2(/(0) + 4)E C 2 - 8] + 
8^2 M S C (E 2 - 1) e/'(0) + 4/(0)E 2 (E 2 - 2) = 0. 

Labels Eigenvalues 

Ef {-2,2,0,0,-^7(07^3^ - 1, vTMVSOxc - l} 

£2 {0,0,0,0,-10 - cos(4u), §0 - cos(4w)| 

E3 {0, 0, -91(7, Qxc, S c ), Si (7, f^xc, 2c), -92(7, ^Xc, 2 C ), 92 (7, Hxc, 2 C )} 

E 4 {0, 0,0,0, -V-97 2 + 9 (7 - 2 )(7 ~ !)2 2 + I87 - 8, V~ 9 7 2 +9(7- 2) (7 - 1)£ 2 + I87 - 8} 

E 5 |0,0, -y/gz{l, X c ,£l X c, 2 C ), ^93(7, ^c.^Xc, 2 C ), -V6s*x c , -y/6x c f (s*)} 

E| {0, 0, Ai(e, E ), A 2 (e, E c ), A 3 (e, S c ), A 4 (e, E c )} 

{0, 0, /J,i(e, 2 C ), At 2 (e, 2 C ), Ms(e, 2 C ), /U 4 (e, 2 C )} 

E 8 {0, 0, -94(7, ^c, Hjc), 94(7, x c, f^Xc), -\/6s*:r c , —V^Xcf (s*)} 

#9 {0, 0, -95(7, ^Xc), 95(7, ^Xc), -96(7, ^Xc), 96(7, ^Xc)} 

Bfo {0. 0- -^^2^*2^97(7, 2 C ), -97(7, 2 C ), -ex/2^/^E c /' («•)} 

E11 {0, 0, -0 (6 - 9E 2 ) + 10S 2 - 4, ^7 (6 - 92 2 ) + 10S 2 - 4, - ^CTv7( 0)E £1 j 

£12 {0, 0, -98(7, 2 C ), 98(7, 2 C ), -99(7, 2 C ), 99(7, 2 C )} 



• Q n is stable for 7 > 1. 

• Qti( s *) is stable for § < 7 < 2, s* > 2, /' (s*) > 0; 
whereas, the possible past attractors are: 

• Qi{s*) is a source for cither < 7 < f , s* < — ^/Jy, /' (s*) < or < 7 < f , s* > y^, /' (s*) > 0; 

• Q5 ( s *) is a source for cither < 7 < |, -\/37 < s* < 0, /' (a*) < or | < 7 < 2, -2 < s* < 0, /' (s*) < or 
0< 7 <|,0<s*< \/37,/' (**) > or | < 7 < 2,0 < s* < 2, f (a*) > 0; 

• Qg and the line Q^[ are unstable for /(0) > 0. 

• is unstable for 7 > 1. 

• Qi2( s *) is source for | < 7 < 2, s* > 2, /' (s*) > 0; 



13 



B. Stability of static solutions 

Having commented the stability issue for expanding (contracting) solutions, let us examine the stability of static 
solutions. As we commented before, the system ([2"2[) - (r26l) admits eighteen classes of (curves of) fixed points corre- 
sponding to static solutions, i.e., having Q = 0, that is with H — 0. Their coordinates in the phase space and their 
existence conditions are displayed in table ITVl 

Now let us comment on the stability of the first order perturbations of (|2ip - ([2"o]) near the critical points showed in 
table [CO 

The line of fixed points E\ is non-hyperbolic. However, the points located at the curve behaves as saddle points 
since its matrix of perturbations admits at least two real eigenvalues of different signs. 

The one-parametric line of fixed points E2 has two real eigenvalues of different signs provided cos(4«) ^ 1. In this 
case, although non-hyperbolic, it behaves as a set of saddle points. For u <G {-| , 4^} all the eigenvalues are zero. In 
this case we need to resort to a numerical elaboration. 

(71(7, S) is always real valued for the allowed values of the phase space variables and the allowed range for the 
free parameters. Then, although non- hyperbolic, the 2D set of fixed points E3 is of saddle type. 

The eigenvalues associated to the one-parametric curve of fixed points E4 are always reals for the allowed values of 
the phase space variables and the allowed range for the free parameters. Thus, although non-hyperbolic, it behaves 
as a saddle point. 

For the allowed values of the phase space variables and the allowed range for the free parameters 33(7, x c , £l\ c , E c ) > 
0. If it is strictly positive, then the fixed points located in the 2D invariant set E$ behaves as saddle points. 

The eigenvalues of the perturbation matrix associated to Eq are {0, 0, Ai(e, S c ), A2(e, E c ), Aa(e, E c ), A4(e, S c )} where 
Ai(e, S c ), A2(e, S c ), As(e, S c ) and A,i(e, E c ) are the roots of the polynomial equation with real coefficients: 

A 4 + A ; V4 - 6£2 e /'(0) + A 2 [/(0) - 4] - 4A^4 - 6E?e/'(0) - 6/(0)£ 2 = 0. (31) 

• For /'(0) 7^ 0, /(0) < 0, equation ([31]) has only two changes of signs in the sequence of its coefficients. Hence, 
using the Descartes's rule, we conclude that for this range of the parameters, there is zero or two positive roots 
of ([3Tj). Substituting A by —A in (|3"Tj) and applying the same rule we have zero or two negative roots of ([31~]) for 
/'(0) 7^ 0, /(0) < 0. If all of them are complex conjugated, we need to resort to numerical investigation. If none 
of them are complex conjugated, then the curve of fixed points Eq consists of saddle points. 

• For e/'(0) > 0, /(0) > 0, equation (f3l"|) has only one change of sign in the sequence of its coefficients. Hence, 
using the Descartes's rule, we conclude that for this range of the parameters, there is one positive root of 
([3"TT) . Substituting A by —A in (f3"Tj) and applying the same rule we have three negative roots of (f5T|) for 
e/'(0) > 0, /(0) > 0. In summary, for e/'(0) > 0, /(0) > 0, at least two eigenvalues of the linear perturbation 
matrix of Eq arc real of different signs. In this case the curve of fixed points Eq consists of saddle points. 

• For e/'(0) < 0, /(0) > 0, equation (|3"Tj) has three changes of sign in the sequence of its coefficients. Hence, 
using the Descartes's rule, we conclude that for this range of the parameters, there are three positive roots 
of (|3"Tj) . Substituting A by —A in ([31]) and applying the same rule we have only one negative root of (|31l) for 
e/'(0) < 0, /(0) > 0. In this case the curve of fixed points Eq consists of saddle points. 

• For f'(0) = /(0) = the non null eigenvalues are —2,2; for /'(0) = 0, /(0) 7^ 0, the non null eigenvalues 

v /- v /iI7(S)S? + (/(0)-4) 2 -/(0)+4 J A /24/(0)El+(/(0)-4)2-/(0)+4 

are i-* — ^= , ±J-2! -j= ; and for / (0) / 0, ./(0) = 0, the non null 

eigenvalues are — 2, 2, -ey/4 - 6£2/'(0). Thus, in both cases consists of saddle points. 

The eigenvalues of the perturbation matrix associated to Ey are {0, 0, /Ui(e, £ c ), ^2(^1 £c), A*3(e, £c)> M4( e > ^c)} where 
/ii(e, E c ), /i2(e, S c ) ) M3( e ; £c) and /(/4(e, S c ) are the roots of the polynomial equation with real coefficients: 

^ + v/2/z 3 S c e/'(0) + fi 2 [2(/(0) + 4)S C 2 - 8] + 8V2 M E C (S 2 - l) e/'(0) + 4/(0)S 2 (S 2 - 2) = 0. (32) 

• For e/'(0) < 0,/(0) < 0,-^ < S c < or e/'(0) > 0,/(0) < 0,0 < S c < or e/'(0) < 0,/(0) < 0,0 < S c < 

or e/'(0) > 0, /(0) < 0, — A= < E c < 0, the equation ([32]) has only two changes of signs in the sequence of 
its coefficients. Hence, using the Descartes's rule, we conclude that for this range of the parameters, there is 
zero or two positive roots of (|32|) . Substituting [i by — fj, in (|32l) and applying the same rule we have zero or 
two negative roots of ([32]) for the same values of the free parameters. If all of them are complex conjugated, 
we need to resort to numerical investigation. If none of them are complex conjugated, then the curve of fixed 
points Ey consists of saddle points. 
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• For /(0) > 0, e/'(0) > 0, < £ c < ^= or /(0) > 0, e/'(0) < 0, ^= < £ c < 0, equation (J32J) has only one change 
of sign in the sequence of its coefficients. Hence, using the Descartes's rule, we conclude that for this range of 
the parameters, there is one positive root of (|32j) . Substituting fi by — fi in (|32j) and applying the same rule we 
have three negative roots of (f3"2"|) for the same values of the parameters. In this case the curve of fixed points 

consists of saddle points. 

• For /(0) > 0, < E c < -j=, e/'(0) < or /(0) > 0, < £ c < 0, e/'(0) > equation (02D has three changes of 
sign in the sequence of its coefficients. Hence, using the Descartes's rule, we conclude that for this range of the 
parameters, there are three positive roots of (|3"2"j) . Substituting [i by —fi in (|3"2"j) and applying the same rule we 
have only one negative root of (|3"2")l for the same values of the parameters. In this case the curve of fixed points 
Ej consists of saddle points. 

• For /'(0) = /(0) = the non null eigenvalues are 2^2^/1 - £j, -2y/2^l - Ef; for /'(0) = 
0,/(0) ^ 0, the non null eigenvalues are ±y / -(/(0) + 4)E 2 - + 4 ) + 16 ) S c " 32E 2 + 16 + 4, and 

±\/-(/(0) + 4 ) S c + v / CT(°)(/(°) + 4 ) + 16 ) S c - 32E 2 + 16 + 4; and for /'(0) ^ 0, / (0) = 0, the non null eigen- 
values are -V^E c e/'(0),2v / 2-2E2, -2^/2 - 2£ 2 . Thus, in both cases £7 consists of saddle points. 

Observe that 54(7, x, fix) is always real valued for the allowed values of the phase space variables and the allowed 
range for the free parameters. Then, although non-hyperbolic, Eg behaves as a set of saddle points. 

When Q\ c 1, 35(7, £l\ c ) is real valued, in this case the fixed points in the line Eg behaves as saddle points. 

For E c ^ 0, 37(7, f^Ac) is real valued, in this case the fixed points in the line E\q behaves as saddle points. 

For E c ^ 0, the linear perturbation matrix evaluated at En has at least two real eigenvalues of different signs, thus, 
the fixed points in the line En behaves as saddle points. 

Observe that the line Eg and the line En contains the special point with coordinates F : (Q = 0,x = 0, VL m = 
0, ft\ = 1,E = 0) in the first case when H,\ c 1 whereas in the second case for E c — s- 0. Taking in both cases 
the proper limits we have that the eigenvalues of the linearization for F are {0, 0, -V67 - 4, ^67 - 4, 0, 0} . Hence, 
for 7 > |, F is of saddle type, whereas for 7 < |, there are two purely imaginary eigenvalues and the rest of the 
eigenvalues are zero. In this case we cannot say anything about its stability from the linearization and we need to 
resort to numerical inspection. 

Observe that gs(j, E) is real valued for the allowed values of the phase space variables and the allowed range for 
the free parameters. Thus, all the points located at the line E12 behaves as saddle points. 

For < 7 < |, Ef s behaves as a saddle point. 

Summarizing, our system admits a large class of static solutions that are proved to be typically of saddle like. This 
kind of solutions arc important as intermediated stages in the evolution of the universe allowing to the transition from 
expanding to contracting models. 



V. PHYSICAL DISCUSSION 



In order to characterize the cosmological solutions associated to the singular points it will be helpful to have the 
important observational parameters in terms of the state variables. The dimensionless scalar field energy density 

parameter fl^ = 3^7, the equation of state parameter = and deceleration parameter q = — ^1 + -jj^j ■ They 

read: 

1 - n A - n m - e 2 2x 2 

a * = Q2 . w * = - 1+ i-n A -n m -£* (33) 

_ PT 2x 2 + (7+l)an + ^A + E 2 -l 

w eff= — = : ^ (34) 

Pt 1 - "a - E^ 

and 

3/ 2»a N ( 2 , 7»m\ , - (. 1V3E 2 

Observe that these parameters, but f2^, blow up as the singular surfaces f^A + Q m + E 2 = 1 and 57a + E 2 = 1 are 
approached. Thus, we need to take the appropriate limits for evaluating at the singular points Qg and Qn- In 
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TABLE VI: Values of the observable parameters ui^, w e ff, and q evaluated at the critical points of the system (|21|l - (|26|l . 
We use the notations a* for an s- value such that f(s*) = 0. 



Pi 


- -< , 




W e ff 


q 


Qf 





arbitrary 


7 


3-1-2 
2 


Qt(s*),Q^(s*) 


i 


1 


1 


2 


Q^(s*) 


s :!: - 


7-1 




2? - 1 


Qs(s*) 


1 


3 V / 


i (s* 2 - 3) 

3 V J 


1 (s* 2 - 2) 

2 V / 







arbitrary 


arbitrary 


2 


Qr 





arbitrary 


arbitrary 


2 


Qt(s*) 


2 / \ 

cos (u) 


1 


1 


2 




1 


-1 


-1 


-1 


Qio 


i-n x 


-1 


-1 


-1 


Qfi 





arbitrary 


arbitrary 


-1 


Q«(0 


4 
s* y 


i 

3 


i 

3 


1 



tabic IVTl we present the values of the observable parameters f2^, w^, u; c ff, and evaluated at the critical points of the 
system (j2"Tj) - (j2l))) that are associated to expanding (contracting) solutions. Static solutions were discussed in section 

EYU 

Observe that, since the variable Q defined in (|T7|) is the Hubble scalar divided by a positive constant, Q > 
corresponds to an expanding universe, while Q < to a contracting one. Furthermore, as usual, for an expanding 
universe q < corresponds to accelerating expansion and q > to decelerating expansion, while for a contracting 
universe q < corresponds to decelerating contraction and q > to accelerating contraction. Lastly, critical points 
with £ = correspond to isotropic universes. 

Now, let us comment on the physical interpretation of the critical points Qi that are associated to expanding 
(contracting) solutions. 

The line of singular points represent a matter-dominated solution (f2 m = 1), as expected they are transient 
stages in the evolution of the universe. 

The singular points Q^(s*) and Q^(s*) are solutions dominated by the kinetic energy of the scalar field. At these 
solutions, the scalar field mimics stiff matter, and tx a -6 where a denotes the scale factor. This solutions represent 
transient states in the evolution of the universe. 

The singular point Q^(s*) represent matter-scalar field scaling solutions that are relevant attractors for either 
< 7 < |,s* < -\/37,/'(s*) < or < 7 < |,s* > yffi, f'(s*) > 0. They satisfy ~ fl m thus they are 
important for solving or alleviating the coincidence problem. The corresponding solutions are accelerating only for 
7 < |. Thus, for standard matter (7 > 1), that is satisfying the usual energy conditions, they cannot represent 
accurately the present universe, since they are decelerating solutions. However, the most interesting feature is that 
the singular point Q4 (s*), which corresponds to a contracting universe (since sgn(Q) = sgn(H) = —1), is a local 
source under the same conditions for which Q~l(s*), which corresponds to an expanding universe, is stable. Therefore, 
in the scenario at hand the transition from contracting to expanding universes is possible, since for the same values 
of the parameters we have a contracting local source and an expanding local attractor. Such a behavior has very 
important physical implications. 

The singular point Q§(s*) is the analogous to P5 in 24]. It represents a scalar- field dominated solution (fl^ = 1) 
that is a relevant attractor for either < 7 < §, -a/37 < s* < 0, /' (s*) < or | < 7 < 2, -2 < s* < 0, /' (s*) < or 
0<7<|,0<s*< \/37;/'(s*) > or | < 7 < 2, < s* < 2, f (s*) > 0. These solutions can represent accurately 
the late-time universe if additionally s* 2 < 2, in which case they represent accelerating solutions. As before, the 
singular point Q~(s*) is a local source under the same conditions for which Q^(s*) is stable. This means that we 
have a large probability to have a transition from contracting to expanding solutions. This is one of the advantages 
of the present model. 

The singular point Q 9 represents an standard 4D cosmological solution whereas Qib represent solutions with 5D- 
corrections. In both cases, from the relationship between y and fl\ follows that this solution is dominated by the 
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potential energy of the scalar field pr = V((f)); that is, it is de Sitter-like solution (w^ = —1). In this case the 
Fricdmann equation can be expressed as 

3H 2 = V (l + |j (35) 

In the early universe, where A <C V, the expansion rate of the universe for the RS model differs from the general 
relativity predictions 

S " /I < 36) 

Contracting de Sitter solutions are associated to the non- hyperbolic fixed point Qg , and it corresponds to the early- 
time universe since it behaves as a local source in the phase-space. 

The singular points Q\ x represent a ID set of singular points, corresponding to isotropic solutions with 5D corrections 
= 1), parameterized by the values of Q = e = ±1 and s c . It can represent the early (late)-time universe for e = +1 
(e = -1) and 7 > 1. 

The singular points Qi 2 (s*) corresponds to a scalar field-dark radiation scaling solution (flu ~ 0^), with an effective 
equation of state parameter of the total matter w e / f = \ corresponding to radiation. These critical points corresponds 
to a relevant late-time attractor for | < 7 < 2, s* > 2, /' (s*) > 0. 

For all the critical point listed above the isotropization has been achieved. The existence of such late-time isotropic 
solutions (e.g., Qt(s*),Q^(s*),Qg , and Qi ), than can attract an initially anisotropic universe, arc of significant 
cosmological interest and have been obtained and discussed in the literature [U, H^- The fact that an isotropic 
solution is accompanied by acceleration, makes it a very good candidate for the description of the observable universe. 

On the other hand, there are solutions that are anisotropic. For example, the solution Qq, and Qy represent 
shear-dominated solution which are always decelerating. Q§ cannot represent accurately the late-time universe, but 
it has a large probability to represent the early-time universe. Qg is of saddle type. For we have similar results 
as for Qq. 

The circles of critical points Qg correspond to scalar field-anisotropic scaling solutions (f2^ ~ E). Both solutions 
represents transient states in the evolution of the universe. When x — > the anisotropic term in the Friedmann 
equation ([5]) dominates the cosmological dynamics (S = ±1). In this limit the corresponding cosmological are in a 
vicinity of the initial singularity. 

Summarizing, the possible late time attractors are: 

• Qt (s*) is stable for cither < 7 < |, s* < -V^Ti /' ( s *) < or < 7 < f , s* > /' ( s *) > 0; 

• Qt(s*) is stable for cither < 7 < f,-\/37 < s* < 0,/'(s*) < or f < 7 < 2, -2 < s* < 0,f (s*) < or 
0<7<§,0<s*< V&yJ' (»*) > or § < 7 < 2,0 < s* < 2,/' (s*) > 0; 

• Qg and the line Qf are stable for /(0) > 0. However, if we include the variable y = jjjt, in the analysis, the 
solutions associated to the line of fixed points Qf are unstable to perturbations along the y-direction. 

• Qi ± is stable for 7 > 1. 

• Qf 2 (s*) is stable for f < 7 < 2, s* > 2, /' (s*) > 0; 
whereas, the possible past attractors are: 

• Qa ( s *) is a source for either < 7 < |, s* < -\/3~7, /' («*) < or < 7 < |, s* > y/5y, f (s*) > 0; 



Q 5 (s*) is a source for either < 7 < |, — s/fry < s* < 0, /' (s*) < or | < 7 < 2, -2 < s* < 0, /' (s*) < or 
0<7<|,0<s*< V5y,f'(s*) > or I < 7 < 2,0 < s* < 2,/' (s*) > 0; 



• Q g and the line Q w are unstable for /(0) > 0. 

• is unstable for 7 > 1. 

• Ql2( s *) is source for | < 7 < 2, s* > 2, /' (s*) > 0; 
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Additionally, our system admits a large class of static solutions that are typically of saddle type. This kind of 
solutions are important at intermediated stages in the evolution of the universe, since they allow to the transition 
from expanding to contracting models and vice versa. In fact, the more interesting result for the scenario at hand is 
that still a large probability for the transition from contracting to expanding universes, since for the same values of 
the parameters we have a contracting local sources and an expanding local attractors. The possibility of re-collapse in 
Bianchi I brane worlds has been investigated for example inJ44| , but in this case for a nonstandard effective equation 
of state of a brane matter p = (7 — l)p with 2 < 7 < 4. In [36[ the authors used an argument similar to Wald's no-hair 
theorem (88j for global anisotropy in the brane world scenarios. There were derived a set of sufficient conditions which 
must be satisfied by the brane matter and bulk metric so that a homogeneous and anisotropic brane asymptotically 
evolves to a de Sitter spacetime in the presence of a positive cosmological constant on the brane. Following that 
reference, the isotropy is reached for any initial condition corresponding to U > (we have obtained the same result 
in (2~H for a quintessence field with positive potential trapped on the brane by without using Wald's arguments). 
When U < the brane may not isotropize and can instead collapse. The existence of re-collapsing solutions have 
been also discussed in [l9|, [2l| for the case of an scalar field with exponential potential on the brane. In the present 
paper we go a step forward by considering potentials beyond the exponential one. 

Finally, a very interesting feature of our scenario is that it allows for bouncing and turnaround, which leads to 
cyclic solutions H§49^. In fact, using a heuristic reasoning we have obtained several conditions for the existence 
of bounces and turnaround of cosmological solutions for a massless scalar field and a scalar field with a potential 
asymptotically constant. For the analysis of more general potentials we submit the reader to the reference [llj where 
the authors discuss conditions for bounces and they provide details of chaotic behavior in the isotropic brane model 
containing an scalar field for a negative dark radiation U < 0. Similarly, in [T|| the results correspond to a perfect 
fluid with equation of state p = (7 — l)p whereas in [98| a scalar field on the brane is considered. Bouncing solutions 
are found to exist both in FRW i?"-gravity_[9J] , as well as in the Bianchi I and Bianchi III i?"-gravity [85| (see also 
[95| ) and in Kantowsky- Sach s i?" -gravity [86|. More interesting, these features can be alternatively obtained in the 
brane cosmology context (96l - ll00| . 

The complete bouncing and cyclic analysis in our framework (that goes beyond the heuristic results presented by 
us in section [TTT|) . as well as the corresponding perturbation investigation, will not be considered in the present work 
and it is left for a future project. Strictly speaking, to complete the bounce and cyclic analysis we have to additionally 
examine the second Friedman equation, that is the equation for H, since one could have the very improbable case 
of H from positive to become exactly zero and then positive again, without obtaining negative values at all. In this 
case we acquire H = and H — simultaneously, that is a universe that stops and starts expanding again without a 
turnaround. 



VI. EXAMPLES OF PHYSICAL POTENTIALS IN THE COSMOLOGICAL LITERATURE 

In this section we want to discuss two physical potentials: the cosh-like potential and the inverted sinh-like potential 
that have been widely investigated in the literature. 

The cosh-like potential has been widely studied for example in [H, [60| - [66| . In [g6[ , and independently by , was 
used a potential of the form V((j>) = Vp [cosh (£0) — 1] to explain the core density problem for disc galaxy halos in the 
ACDM model (see also [48[ and [iij], HH). The asymptotic properties of a cosmological model with a scalar field with 
cosh- like potential have been investigate in general relativistic framework by [46| and for RS2-FRW branes in [2^| . 

The inverted sinh-like potential has been widely investigated in the literature for example in [4^, 0, [6^, [(37| ■ The 
asymptotic properties of a cosmological model with a scalar field with such a potential have been investigated in the 
context of FRW brane in [H]. 

A. Potential V{<j>) = V [cosh - 1] 

Now let us consider the potential 

F(0) = Vb[cosh(^)-l] (37) 
The /-diviser corresponding to this potential is given by 

f{s) = -\{s 2 -e) (38) 

which has its zero is 

s*=±i /V) = T£ (39) 
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(b) 



FIG. 1: (a) Projection of some orbits of the system (|2HI - (|26| I for the potential V((f>) = Vo [cosh (£0) — 1] with £ = 1/2 in the 
planes x — Q.\ and Vl a — Q,\ for 7 = 1. This numerical elaboration shows that Qg is the local attractor. 

(b) Projection of some orbits in the space E, Q, Q. m . We choose the initial states near static solutions. This numerical elaboration 
shows the transition from the contracting de Sitter solution Qg (past attractor) to the expanding de Sitter solution Qg (future 
attractor) . 



The only possible expanding late time attractor of the model with potential (|37[) is the Sitter attractor Qg 

In figure [1] are showed some numerical integrations for the system (f2"Tj) - (j2"(3)) for the function (f55|) with £ = 1/2. This 
numerical elaboration shows that the expanding de Sitter solution Qg is the future attractor whereas the contracting 
de Sitter solution Qg is the past attractor. It is illustrated also the transition from contracting to expanding solutions. 
In figure [T](b) we choose initial conditions near the static solutions presented in table HVl This kind of solutions allow 
for the transition form contracting deceleration to expanding accelerated solutions. 

B. Potential V((f>) = V Biah~ a (j34>) 

Let us consider the potential 

V{4>) = VQsmh- a {pcj>) (40) 
The /-diviser corresponding to this potential is given by 

f(s) = --a(3 2 (41) 



19 



The zeros of this f(s) function are 

s*=±a(3 /'(s*) = ±2/3 (42) 
The sufficient conditions for the existence of late-time attractors are fulfilled easily 

• The scalar field-matter scaling solution Qf(a(3) is a late-time attractor provided 0<7< < 0,a > — -^g j 
or < 7 < §,a > ^F,/3 > 0. 

• The scalar field-matter scaling solution Q^(— a/3) is a late-time attractor provided < 7 < |,q < -^|p,/3 > 
or < 7 < f,a > -^f-,P < 0. See figured 

• The solution dominated by scalar field Q~^(a(3) is a late-time attractor provided < 7 < |, /3 < 0, < a < — 

or I < 7 < 2,/3 < 0,0 < a < -§ or < 7 < |,0 < a < ^-,(3 > or § < 7 < 2, < a < f , /3 > 0; 

• The solution dominated by scalar field Q^(— a/3) is a late-time attractor provided < 7 < |,/3>0, 0<a< 

or|<7<2,/3>0,0<a<|or0<7<|,0<a< -^,/3 < or | < 7 < 2, < a < -|,/3 < 0; see 
figures [3] 

• The de Sitter solution Qt[ is the late time attractor provided /(0) = —a/3 2 > 0. See figure SJ 

• Qi ± is stable for 7 > 1. 

• The scalar field- dark radiation scaling solution Qf 2 (a(3) is a late-time attractor for | < 7 < 2, a/3 > 2, /3 > 0. 

• The scalar field- dark radiation scaling solution Qf 2 (—a(3) is a late-time attractor for | < 7 < 2, a/3 < —2, /3 < 0. 
The past attractors of the system are as follows. 

• The scalar field- matter scaling solution Q^(af3) is the past attractor provided 0<7< |, /3 < 0, a > —jp-, or 
< 7 < |,a > > 0. 

• The scalar field-matter scaling solution Q^(— a/3) is the past attractor provided < 7 < |,a < y -^-,f3 > or 
< 7 < |,a > -^,0< 0. See figure [2] (b). 

• The solution dominated by scalar field (a/3) is the past attractor provided 0<7< |,/3<0,0<a< —jp- 
or|<7<2,/3<0,0<a<-|or0<7<|,0<a< ^,(3 > 0or|< 7 <2,0<a<|,^>0; 

• The solution dominated by scalar field Q$(— a/3) is the past attractor provided < 7 < |,/3 > 0, 0<a< 

or|<7<2,/3>0,0<a<|or0<7<|,0<a< -^,/3 <0or|<7<2,0<a< -|,/3 < 0; see the 
figure 1(b). 

• The de Sitter solution Q$ is the past attractor provided f(0) = —a/3 2 > 0. See figure 0] (b) 

• is unstable for 7 > 1. 

• The scalar field- dark radiation scaling solution Qj" 2 (a/3) is a past attractor for | < 7 < 2, a/3 > 2,/3 > 0. 

• The scalar field- dark radiation scaling solution Qi 2 (— a/3) is the past attractor for | < 7 < 2, a/3 < —2, /3 < 0. 

In figure [2] we present some orbits in the phase space of (|2T |) -([26 |) for the potential V(4>) = Vq s'mh~ a (f3(f>) with 
a = 1/2, P = 4 and 7 = 1. In figure [2] (a) we set Q = +1 for the numerical simulation. This numerical elaboration 
reveals that Q\ (1/8) is the future attractor. In figure [5] (b) we present a projection of some orbits in the subspace 
(S, Q, O m ). We choose the initial states near static solutions. This numerical elaboration illustrates the transition for 
the contracting scalar-field matter scaling solution Qj(l/8) (past attractor) to the expanding one Q~£(l/8) (future 
attractor). As can be seen from the figures, the static solutions play an important role in the transition form contracting 
to expanding solutions and vice versa. In figure [5] (c) we show the transition from expansion to contraction and then 
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(c) 

FIG. 2: Some orbits in the phase space of ((2T]) - (f26l) for the potential V((j>) = Vo sinh~ a (/3<j>) with a = 1/2, f) = 4 and 7 = 1. (a) 
Projection in the planes x — Q\ and x — Q m . We set Q — +1 in the numerical simulation. This numerical elaboration shows 
that Qj(l/8) is the future attractor. 

(b) Projection in the subspace (£,Q,fi m ). We choose the initial states near static solutions. This numerical elaboration 
illustrate the transition for the contracting scalar-field matter scaling solution Qj(l/8) (past attractor) to the expanding one 
Q4 (1/8) (future attractor). 

(c) Transition from expansion (Q > 0) to contraction (Q < 0) and then back to expansion (Q > 0), that is a cosmological 
turnaround and a bounce. The left panel corresponds to the projection of an orbit in the plane x-Q, while the right panel 
shows Q = 10Q vs r, in order to illustrate that Q change its sign twice during the evolution. We choose the initial condition 
Q(0) = 0.9, x(0) = 0.6, Q m (0) = 0.25, fi A (0) = 0, E(0) = 0, s(0) = 2.0. 
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FIG. 3: Some orbits of the phase space of (gTJl-® for the potential V{4>) = Vq sinh" a (/30) with a = /3 = 1/2 and 7 = 1. 
(a) Projection in the planes x — £l\ and Sl a — Q\ (We set Q = +1 in the numerical simulation), (b) Some orbits in the 
projection (E, Q, Q m )- We choose initial states near the static solutions. This numerical elaboration shows the transition from 
the contracting scalar field dominated solution Q^(l/4) to the expanding one Qj|"(l/4). 



back to expansion, that is a cosmological turnaround and a cosmological bounce, for the solutions of ([2"T]) - ([2"6|) . for 
the potential V(<f>) = Vq sinh _Q {13(f) with a = 1/2,(3 = 4 and 7 = 1. The left panel corresponds to the projection 
of an orbit in the plane x-Q, while the right panel shows Q = 10Q vs t in order to illustrate that Q change its sign 
twice during the evolution. We choose the initial condition Q(0) = 0.9, x(0) = 0.6,f2 m (0) = 0.25,J1a(0) = 0, E(0) = 
0,s(0) = 2.0. 

In figure [3] are presented some orbits in the phase space of (|2T1) - ([2^1) for the potential V((f>) = Vq s'mh~ a ((3(f>) with 
a = j3 = 1/2 and 7=1. In figure [3] (a) we set Q = +1 in the numerical simulation and we present the projection 
in the planes x — Q\ and — Q\. In this case the future attractor is the isotropic solution dominated by scalar 
field Q^{\/4), also this solution corresponds to an accelerated expansion rate since s* = 1/4 < \/2. In figure [3] (b) 
are presented some orbits in the projection (E,<3,f2, rl ). This numerical elaboration shows the transition from the 
contracting scalar field dominated solution (1/4) to the expanding one <5^"(l/4). We choose the initial states near 
the static solutions. 

In figures |4] are showed some numerical integrations for the system (|2~Tj) - (f26| for the function (|4T|) with a = 
— 1/2,/? = —1/4. In this case the local attractor is the solution dominated by the potential energy of the scalar field 
Qg (Sitter attractor). The past attractor is Qg . This numerical elaboration illustrate the transition from contracting 
to expanding solutions and vice versa that is not allowed for Bianchi I branes with positive Dark Radiation term 
{U > 0). 
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FIG. 4: Some orbits of the phase space of $T$-$$ for the potential V{4>) = V Q sinh _O! (/30) with a = -1/2,(8 = -1/4 and 
7 = 1. (a) Projection in the planes x — Q\ and E — fl m . We set Q = +1 in the numerical simulation. 

(b) Projection in the subspace (E, Q, Q m ), Q 7^ +1. We choose the initial states near static solutions. This numerical elaboration 
illustrate the transition for the contracting de Sitter solution Qg (past attractor) to the expanding de Sitter solution Qg (future 
attractor) . 



VII. CONCLUSIONS 

In this paper we have presented a full phase space analysis of a model consisting of a quintessence field with 
arbitrary potential and a perfect fluid trapped in a Randall-Sundrum's Brancworld of type 2. We have considered 
a homogeneous but anisotropic Bianchi I (BI) brane geometry. Moreover, we have also included the effect of the 
projection of the five-dimensional Weyl tensor onto the three-brane in the form of a negative Dark Radiation term. 

We have discussed a general method for the treatment of the potential, called "Method of /-devisers" that allows 
investigating arbitrary potentials in a phase space. This method is a modification of the method that has been applied 
to isotropic (FRW) scenarios, and that has been generalized to several cosmological contexts. The "Method of /- 
devisers" has the significant advantage, that one can first perform the analysis for arbitrary potentials and then just 
substitute the desired forms, instead of repeating the whole procedure for every distinct potential. More importantly 
is that the method does not depend on the cosmological scenario, since it is constructed on the scalar field and its self- 
interacting potential, and can be generalized to several scalar fields. The disadvantages are that such a generalization 
is only possible if in the model do not appear functions that contains mixed terms of several of the scalar fields, or if 
the model contains more than one "arbitrary" function, for example non-minimally coupled scalar field cosmologies. 
In such a case, the more convenient approach is to consider the scalar field itself as a dynamical variable. 

By means of the "Method of /-devisers" as well as some tools from the Theory of Dynamical Systems, we have 
obtained general conditions which have to be satisfied by the potential (encoded in the mathematical properties of the 
key /(s)-functions) in order to obtain the stability conditions of standard 4D and non-standard 5D de Sitter solutions. 

We have presented general conditions under the potential for the stability of standard 4D de Sitter solutions. We 
proved that the de Sitter solutions with 5D corrections are stable against the perturbations introduced here, but are 
unstable against small perturbations of V/(3H 2 ). This fact does not conflict with our previous results for BI branes 
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with positive dark radiation term. Also, we presented the stability conditions for both scalar field-matter scaling 
solutions, for scalar field-dark radiation scaling solutions and scalar field-dominated solutions. The shear-dominated 
solutions are always unstable (contracting shear-dominated solutions are of saddle type). 

A crucial difference with respect our previous works is that for hi < the traditionally ever-expanding models 
could potentially re-collapse due to the of this dark radiation. Additionally, our system admits a large class of static 
solutions that are of saddle type. This kind of solutions are important at intermediate stages in the evolution of 
the universe, since they allows to the transition from contracting to expanding models and vice versa. Thus, a new 
feature of this scenario is the existence of a bounce and a turnaround, which leads to cyclic behavior. This features 
are not allowed in Bianchi I branes with positive dark radiation term. This is an important cosmological result in our 
scenario, that is a consequence of the negativeness of the dark radiation. 

Finally, in order to be more transparent, we have illustrated our main results for the specific potentials V oc 
sinh _a (/3(f)) and V oc [cosh — 1] which have simple /-devisers. 
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Appendix A: De Sitter solutions. 

To investigate the de Sitter solution Qg , we introduce the local coordinates: 

{x, fl x ,Q - 1, fi m , E, s] = e {x, fl\,Q, fl m , &*,s} + O(e) 2 (Al) 
where e is a small constant eCl. Then, we obtain the approximated system 

Q' = -4Q,x' = y|s-3f,Q' A = 0, 

Q' m = -3 7 Q ro , Cl' a = -3tl a , s' = -V6f(0)x. (A2) 
The system (|A2[) admits the exact solution passing by (Qo, xq, fl\o, fl m o, ^o-O; so) at t = given by 
Q(r) = Q e- 4r , 

x(r) = \x e-^ (e- + l) - ^ = = 9 + ^ , 

Clx{r) = Oao, 

m (r) = Q m0 e- 3 ^,n a (T) = £W~ 3r , 

* io , x 3e~ 3r / 2 (2s„ - V6x ) sinh(/3r) j~2 , /0 
J(r) = s e-* T/2 cosh(/?r) + ^ ~ + V ' ^HM (A3) 



where /3 = ^9-12/(0). 

For the choice /3 2 < |, i.e., /(0) > 0, all the perturbations Q, x, Cl m , Cl a , s shrink to zero as r — > +oo, and fl\ 
converges to a constant value. For (3 = ±|, i.e., for /(0) = 0, x — > and s — > sq, and the other perturbations shrink 
to zero as r — > +oo. Combining the above arguments we obtain that for /(0) > 0, Qt is stable, but not asymptotically 

32 ^ 9 



stable. For /3 > |, i.e., /(0) < 0, all the perturbation values, but x and s, which diverges, shrink to zero as r — > +oo. 
+ 

9 



Thus, Q+ is a saddle. 
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For analizing the curve of singular points Qio we consider an arbitrary value f2 A = u c , < u c < 1, and introduce 
the local coordinates 

{x, n x - Uc, Q - 1, On, E, s} - e {£, Oa, 4 n m , fi ff) s} + 0(e) 2 (A4) 

where e is a small constant e <C 1. 

Then, we obtain the approximated system 

Q' = -4Q, £' = -3x - y|s(u c - 1), Cl' x = -3ju c £l m , 

Q' m = -3 7 f2 ro , n; = -30 CT , s' = -VEf(0)x. (A5) 
The system (| A5|) admits the exact solution passing by (Qo, xq, ^ACb ^mOi &<ro, so) at r = given by 

Q(t) = Qoe ,Jc(t) = — he aT/z a;o cosher), 

fi A (r) = u c Q m0 (e- 3 ^ T -l) + O A0 , 

fim(r) = a„ e- 37T ,^(r) = fW~ 3r , 

_ 3e -H^) (e^-1) (23o( Ue -l) + V6So) 1 , ^ lU 

[) ~ 4(u c -l)0 2 S ° 1 + j + 



£0^-^+3) (e T * - 1) 
2\/6K - 1) 



(A6) 



where d = ^/l2/(0)(u c — 1) + 9. It is easy to see that for /(0) > the solution is stable, but not asymptotically 
stable. For /(0) < is of saddle type. 

It is worthy to mention that the above results could be proved by noticing that Q^ -as ID set- is normaly hyperbolic 
since the eigen-direction associated with the null eigenvalue, (0, 0, 1, 0, 0, 0) T , that is, the f2 A -axis, is tangent to the 
set. Thus, the stability issue can be resolved by analyzing the signs of the non-null eigenvalues. However, let us 
comment that if we include the variable y = gjjr, in the analysis, the solution is unstable to perturbations along the 
^-direction. 

Appendix B: Asymptotic analysis on the singular surface S7a + £ 2 = 1. 

In this section we analyze the asymptotic structures of the system (f2"Tj) - (|2l)]) near the critical points at the singular 
surface fl x + E 2 = 1. 

1. Shear-dominated solution. 

Observe that the one of the singular points at the surface f2 A + Y? = 1 that represents a shear-dominated solution 
is Qg. To investigate its stability we introduce the local coordinates: 

{x, n x ,Q - 1, n m , s + 1, s - s c } = e {x, n Xl q, n m ,n a , s} + o(e) 2 (bi) 

where e is a constant satisfying e << 1 and s c is an arbitrary real value for s. Taylor expanding the resulting system 
and truncating the second order terms we obtain the approximated system 

q' = 2Q,x' = - y|s c (-2Q CT + il x + n m ) , 

«A = 60 A ( 7 °"' + 1 ) , Cl' m = 3(2 - 7 )f2 m , 
y n x — 2ii a j 

^ = ?^Sr + ^ §> = ~^fM- ( B2 ) 
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Let us define the rates: 

{r m ,r a } = -^{£l m ,n a } (B3) 
"A 

Let us assume 7 > 0. With the above assumptions we can define the new time variable N = 7 In a which preserves 
the time arrow. Then we deduce the differential equations 

6r 2 

3r m (B4) 



m 2r a - 1 
and 

, 3r rr 

r > = — 

where the comma denotes derivatives with respect to N. 
The system admits the general solution: 



1 _ 2ci - 1 + tanh (i (3N - c 2 )) 

-2ci cosh (37V -c 2 )-2ci' r,T ~ 4^ 

1 _ 2d - 1 - tanh (| (3iV - c 2 )) 



(B5) 



(B6) 



(B7) 



-2ci cosh (3JV - c 2 ) - 2ci 4ci 

where Ci 7^ and c 2 are arbitrary constants. In the general case we have that r m — > and tends to a constant as 
N -> -00. 

Thus, as r — > —00 the system (|B2I) has the asymptotic structure: 

Q' = 2Q, f = -w|s c (f2 A - 2fi a ), 0' A = 60 A , 

o;„ = 3(2 - 7 )n m , = er> CT , s' = - V6xf( Sc ) (B8) 

The system (|B8[) admits the exact solution passing by (Qo, £o 5 ^A(h ^mOj ^o-Oj ^0) at r ~ given by 

Q(r) - Q e 2T , sg(t) = £ - — , fl x (r) = tt A0 e 6r , 

Cl m (r) = a„oe 3(2 " 7)T ,^(r) = Cl aQ e 6T , 

S(t) = Y^s c (-6r + e 6T - l) /(s c )(f2 A0 - 2ft ff0 ) - Vftr£o/(«c) + s . (B9) 

Thus, the energy density perturbations goes to zero, and s diverge as r — > —00. For that reason, Q§ cannot be the 
future attractor of the system, thus Qq has a large probability to be a past attractor. Using the same approach for 
Qg we obtain that as r — > +00, the energy density perturbations goes to zero, but s diverge. Thus, although Q§ 
have a large probability to be a late time attractor, actually it is of saddle type. For we have similar results as 
for Q 6 . Since the procedure is essentially the same as for Qg , we omit the details. 



2. Solution with 5D-corrections. 

Another singular points at the surface £l\ + S 2 = 1 are the isotropic solutions with Q\ = 1 (critical points Qfi). 
For the stability analysis of , we introduce the local coordinates: 

{x, Q\ — 1,Q — l,n m , S, s - s c } = e {x, Cl x , Q, fi m , f2 CT) s} + 0(e) 2 (BIO) 

where e is a constant satisfying e << 1 and Q c , s c are arbitrary real values of Q, s respectively. Then, we obtain the 
approximated system 

Q' = -2Q(3 7 r + 2), x' = - (-6z( 7 r + 1) - V6s c (0 A + O m )) , 
Cl' x = -3 7 rf2 A , f2' m = -3 7 (2r + l)Q m , 

tt' a = -36 CT ( 7 r + 1), s' = -V6£/(s c ). (Bll) 
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where we have defined the rate r = The evolution equation for r is given by 

r' = -3jr(r + 1). (B12) 

The equation (|B12j) has two singular points r = with eigenvalue ^jp\ r =o = ~37 and r = — 1 with eigenvalue 
^p| r =-i = 37. Then, assuming 7 > we have that r — > —1 as r — > —00. This argument allow us to prove that as as 
r — > —00, the system (|B11[) has the asymptotic structure 

Q' = 2(3 7 - 2)Q, x' = 35(7 - 1), n' A = 3 7 fi A , 

fi^ = 3 7 a„, fi^ = 3^(7 - 1), s' = -V6xf(s c ). (B13) 

The system (|B13|1 admits the exact solution passing by (Qq, Xq, Cl\ Q , tl m0 ^ Cl^Q, sq) at r = given by 

Q(r) = Q e 2 ^- 2)T ,x( T ) = x e 3 ^ T , Cl x (r) = Cl X0 e 3 ^, 
O m (r) = Cl m0 e^ T ,n a (r) = jW 3(7 " 1)T , 

J\x f{s c ) {e^-^ 1) 
«(r) = h ^ (B14) 

Hence, for 7 > 1, the perturbations Q, x, f2 m , shrink to zero as r — > —00. In that limit s tends to a finite value. 
In this case, the solution under investigation has a 4D unstable manifold. The numerical simulations suggest that the 
singular point at the surface tt\ + E 2 = 1 with £l\ = 1 is a local source. On the other hand, for 7 < 1, the solution 
behave as a saddle point, since the perturbation values x and Q, a do not tend to zero as r — > —00. Using the same 
approach for we obtain that for 7 > 1, is stable as r — > +00, whereas for 7 < 1 it is of saddle type. These 
results for and for are in agreement with the analogous results for m+ and for ?n_ in the reference (T^ |. 
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